print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(G9a-/- vs. control)")+
theme_classic2(base_size = 20)
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.5))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(G9a-/- vs. control)")+
theme_classic2(base_size = 20)
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(G9a-/- vs. control)")+
theme_classic2(base_size = 20)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(G9a-/- / control)")+
theme_classic2(base_size = 20)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(ratio between G9a-/- & control flies)")+
theme_classic2(base_size = 20)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(ratio between G9a-/- & control flies)")+
theme_classic2(base_size = 16)
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("Controls_riskratio.csv")
hazfn <- read.csv("Controls_riskratio.csv")
str(hazfn)
summary(hazfn)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(ratio between G9a-/- & control flies)")+
theme_classic2(base_size = 16)
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("Controls_riskratio.csv")
hazfn <- read.csv("Controls_riskratio.csv")
str(hazfn)
summary(hazfn)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(ratio between outross & y/w flies)")+
theme_classic2(base_size = 16)
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("Risk ratio.csv")
hazfn <- read.csv("Risk ratio.csv")
str(hazfn)
summary(hazfn)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(ration betweem Jak/Stat mutants & y/w)")+
theme_classic2(base_size = 20)
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("akStat_riskratio.csv")
hazfn <- read.csv("JakStat_riskratio.csv")
str(hazfn)
summary(hazfn)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(ration betweem Jak/Stat mutants & y/w)")+
theme_classic2(base_size = 20)
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("JakStat_riskratio.csv")
hazfn <- read.csv("JakStat_riskratio.csv")
str(hazfn)
summary(hazfn)
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(ration betweem Jak/Stat mutants & y/w)")+
theme_classic2(base_size = 20)
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("JakStat_riskratio.csv")
hazfn <- read.csv("JakStat_riskratio.csv")
str(hazfn)
summary(hazfn)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(Jak/Stat mutants vs. y/w)")+
theme_classic2(base_size = 20)
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("JakStat_riskratio.csv")
hazfn <- read.csv("JakStat_riskratio.csv")
str(hazfn)
summary(hazfn)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(Jak/Stat mutants vs. y/w)")+
theme_classic2()
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("JakStat_riskratio.csv")
hazfn <- read.csv("JakStat_riskratio.csv")
str(hazfn)
summary(hazfn)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(Jak/Stat mutants vs. y/w)")+
theme_test()
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("JakStat_riskratio.csv")
hazfn <- read.csv("JakStat_riskratio.csv")
str(hazfn)
summary(hazfn)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(Jak/Stat mutants vs. y/w)")+
theme_classic()
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("JakStat_riskratio.csv")
hazfn <- read.csv("JakStat_riskratio.csv")
str(hazfn)
summary(hazfn)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(Jak/Stat mutants vs. y/w)")+
theme_classic(base_size = 18)
setwd("~/Dropbox/Mickael/Analysis/R")
read.csv("JakStat_riskratio.csv")
hazfn <- read.csv("JakStat_riskratio.csv")
str(hazfn)
summary(hazfn)
## boxplot
library(ggplot2)
ggplot(data = hazfn) +
aes(x = flyline, fill = sex, color = sex, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "Flyline",
y = "Risk ratio (Jak/Stat mutants vs. y/w)") +
theme_test()
p<- ggplot(hazfn, aes(x=flyline, y=estimate, fill=sex)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9))
print(p)
# Finished bar plot
p+labs(x="Flyline", y = "Risk ratio
(ratio between Jak/Stat mutants & y/w)")+
theme_classic(base_size = 18)
setwd("~/Dropbox/Thesis/Chapter-3 & 4/Gut-epithelia & Jak:stat/Project_Jak:Stat-pathway & tolerance/Analysis/R/csv files")
read.csv("Linear_slope.csv")
slope <- read.csv("Linear_slope.csv")
str(slope)
summary(slope)
## boxplot
library(ggplot2)
ggplot(data = slope) +
aes(x = line, fill = line, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "FLy line",
y = "Linear slope")  +
theme_test(base_size = 18)+
facet_grid(cols = vars(sex))
p<- ggplot(slope, aes(x=line, y=estimate, fill=line)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9)) + facet_grid(cols = vars(sex))
setwd("~/Dropbox/Thesis/Chapter-4/Project_Jak:Stat-pathway & tolerance/Analysis/R/csv files")
setwd("~/Dropbox/Thesis/Chapter-4/Project_Jak:Stat-pathway & tolerance/Analysis/R/csv files")
read.csv("Linear_slope.csv")
slope <- read.csv("Linear_slope.csv")
str(slope)
summary(slope)
## boxplot
library(ggplot2)
ggplot(data = slope) +
aes(x = line, fill = line, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "FLy line",
y = "Linear slope")  +
theme_test(base_size = 18)+
facet_grid(cols = vars(sex))
p<- ggplot(slope, aes(x=line, y=estimate, fill=line)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9)) + facet_grid(cols = vars(sex))
print(p)
setwd("~/Dropbox/Thesis/Chapter-4/Project_Jak:Stat-pathway & tolerance/Analysis/R/csv files")
read.csv("Linear_slope.csv")
slope <- read.csv("Linear_slope.csv")
str(slope)
summary(slope)
## boxplot
library(ggplot2)
ggplot(data = slope) +
aes(x = line, fill = line, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "FLy line",
y = "Linear slope")  +
theme_test(base_size = 18)+
facet_grid(cols = vars(sex))
p<- ggplot(slope, aes(x=line, y=estimate, fill=line)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9)) + facet_grid(cols = vars(sex))
print(p)
# Finished bar plot
p+labs(x="Fly line", y = "Tolerance Slope
(estimates from linear reaction norm)")+
scale_fill_manual(values = c("gray27","indianred3","slateblue3","goldenrod","plum3"))+
theme_pubr(base_size = 18)+
ylim(20, -55)
setwd("~/Dropbox/Thesis/Chapter-4/Project_Jak:Stat-pathway & tolerance/Analysis/R/csv files")
read.csv("Linear_slope.csv")
slope <- read.csv("Linear_slope.csv")
str(slope)
summary(slope)
## boxplot
library(ggplot2)
ggplot(data = slope) +
aes(x = line, fill = line, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "FLy line",
y = "Linear slope")  +
theme_test(base_size = 18)+
facet_grid(cols = vars(sex))
p<- ggplot(slope, aes(x=line, y=estimate, fill=line)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9)) + facet_grid(cols = vars(sex))
print(p)
# Finished bar plot
p+labs(x="Fly line", y = "Tolerance Slope
(estimates from linear reaction norm)")+
scale_fill_manual(values = c("gray27","indianred3","slateblue3","goldenrod","plum3"))+
theme_classic(base_size = 18)+
ylim(20, -55)
setwd("~/Dropbox/Thesis/Chapter-4/Project_Jak:Stat-pathway & tolerance/Analysis/R/csv files")
read.csv("Linear_slope.csv")
slope <- read.csv("Linear_slope.csv")
str(slope)
summary(slope)
## boxplot
library(ggplot2)
ggplot(data = slope) +
aes(x = line, fill = line, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "FLy line",
y = "Linear slope")  +
theme_test(base_size = 18)+
facet_grid(cols = vars(sex))
p<- ggplot(slope, aes(x=line, y=estimate, fill=line)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9)) + facet_grid(cols = vars(sex))
print(p)
# Finished bar plot
p+labs(x="Fly line", y = "Tolerance Slope
(estimates from linear reaction norm)")+
scale_fill_manual(values = c("gray27","indianred3","slateblue3","goldenrod","plum3"))+
theme_classic(base_size = 18)
setwd("~/Dropbox/Thesis/Chapter-4/Project_Jak:Stat-pathway & tolerance/Analysis/R/csv files")
read.csv("Linear_slope.csv")
slope <- read.csv("Linear_slope.csv")
str(slope)
summary(slope)
## boxplot
library(ggplot2)
ggplot(data = slope) +
aes(x = line, fill = line, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "FLy line",
y = "Linear slope")  +
theme_test(base_size = 18)+
facet_grid(cols = vars(sex))
p<- ggplot(slope, aes(x=line, y=estimate, fill=line)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9)) + facet_grid(cols = vars(sex))
print(p)
# Finished bar plot
p+labs(x="Fly line", y = "Tolerance Slope
(estimates from linear reaction norm)")+
scale_fill_manual(values = c("gray27","indianred3","slateblue3","goldenrod","plum3"))+
theme_classic(base_size = 20)
setwd("~/Dropbox/Thesis/Chapter-4/Project_Jak:Stat-pathway & tolerance/Analysis/R/csv files")
read.csv("Linear_slope.csv")
slope <- read.csv("Linear_slope.csv")
str(slope)
summary(slope)
## boxplot
library(ggplot2)
ggplot(data = slope) +
aes(x = line, fill = line, weight = estimate) +
geom_bar(position = "dodge") +
labs(x = "FLy line",
y = "Linear slope")  +
theme_test(base_size = 18)+
facet_grid(cols = vars(sex))
p<- ggplot(slope, aes(x=line, y=estimate, fill=line)) +
geom_bar(stat="identity", color="black",
position=position_dodge()) +
geom_errorbar(aes(ymin=estimate-std_err, ymax=estimate+std_err), width=.2,
position=position_dodge(.9)) + facet_grid(cols = vars(sex))
print(p)
# Finished bar plot
p+labs(x="Fly line", y = "Tolerance Slope
(estimates from linear reaction norm)")+
scale_fill_manual(values = c("gray27","indianred3","slateblue3","goldenrod","plum3"))+
theme_classic(base_size = 22)
